LER ES Be os 
; t ts — : 


| March 1982 
; - 
i. co 
ei ti 
1° ae 
. - prepared for: : Goes 
: = : “Waval Postgraduate Schoo! Ace ma a 
ie eee ae | pMonterey, california 93946 “yomait bby 

; on Aa ele "ime. oe La wba iP : 

: Te Mi eae 


ADAI11I15108 


oJ) 


NPS67--82~002CR 


NAVAL POSTGRADUATE SCHOGL 


Monterey, California 


CONTRACTOR REPORT 


-PERFORMARCE OF SOLID FUEL RANJET 
GUIDED PROJECTILE FOR USN 
B'754 GUN SYSTEM 


“ODED AMICHAI 


a > eae So ere een Oy eee ee Se Mr er eee + o Soe: 
: (i 
Ne 
we 
| NAVAL POSTGRADUATE SCHOOL 
Monterey, California 
{ 
Rear Adwiral J, J. Ekelund . D. A. Schrady : 
Superintendent Acting Provost 


The work reported herein was carried out for the Naval Postgraduate 
School by Oded Anichat under Contract Humber 00228-81-~C-H231, The work 
presented in this report is in support of solid fuel ranjet research 
and the exploration of Navy applications for Advanced Indfrect Fire : » 
Support, AIFS. Both projects are funded by the Defense Advanced Research . 
Projects Agency and are under the cognizance of Professor A. €. Fuhs. : 


Reproduction of all or part of this report is author {zed. 


i" 
fy . 
3 \4 _ This report was prepared by: - 
Yj tA : . ia f 
oe Contracted Research ‘Associate 
a > Reviewed by: a a . - : 
eS S BSD Sagat Fy ee r 
i Sa iy aoe , 
el i (tlie Zz CD oreo a! tthe } i 
aa Brien E. % . a ik eld. ns ot 
- Distinguishsd Professor == “Acting Giairaan, Defa. trent of Aeronautics — 
j 


“. : , i ¥, Mw. Tolles 
“Daan of Research 


t 


. ose 


ee 
ie. ce ee en ee ee ee ene 
F =e : es 


h al 
if 


! DISCLAIMER NOTICE 


, THIS DOCUMENT 1S BEST QUALITY 
bas PRACTICABLE. THE COPY FURNISHED 
: TO DTIC CONTAINED A SIGRIFICANT 


uy NUMBER OF PAGES WHICH DO NOT 
: REPRODUCE LEGIBLY. | 


ce a a ee 


— 


BO Seg gree bail > i ae ie 5 See Ty - . 
iy 
i 
i 


SUNT EAS S LEED 


SECUR TY CLASSIFICATION OF THIS PAGE (When Date E.tered) 


RMS REPORT DOCUMENTATION PAGE [pee AD IRRRUC TIONS 


REPOR? “Sen 7. GOVT ACCESSION NQ{ 3. RECIPIENT S CATALOG NUMBUM 
NPS67-62-002CR h-. Aro. 


& TITLE Cand Subticte) 


3. TYPE oF REPORT & PERIOD CCVERED 
PERFORMANCE OF SOLID FUEL RAMJET GUIDED | ‘ 
PROJECTILE FOR USM 5"/54 SUN SYSTEM RAACT aE REROTE 

[© Pearcemia O8G. REPORT NUMBER 


S. TONTAACT OR GRANT NUGRERCE) 
NON228-81-C-H231 


t 


rape eee 


em 0 eetemmeemamer 
AUTHOR! a) 


“(ded Amichai 


7 


en certo een as mee er 
10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMGERS 


DARPA ORDER 4035 
TELEMENT: 62702E 


9 PERFORMING ONGANILE TIOnw NAME AND ADUNESS 


Navai Postcraduate School 
Monterey, California 93940 


emer renlere v 


arenes, 
12. REPORT OATE 


March 1982 


1a, NUMOER.OF F AGES 
saber cy) 


VS. COBTROLLING GT FICE NAME AND ADDRESS 


Defense Advanced Research Projects Agency 
1400 Wilson Boulevard 
Arlington, Virginia 22209 


% MONITORING AGENCY THAME & AOODRESSI(IE different fram Controlling Oftice) 


a NTI HA © ST me 
1S. SECURITY CLASS, (of this report) 


UNCLASSTFIED 


Se. OUCLASSIFIGATION/ DOWNGRADING 
SCrEOULE 


(6. DISTRIBUTION STATEMENT fol thta Report) 


—. 


ee ne 


Approved for public release, distributiesiinlimited, 


a 
17, BISTRIGUTION STATEMENT (at the abetract entacod in Block 10, (f different from Report) 


SUPPLEMENTARY NO-ES 


ee tte nar ne en RR 
19. KEY WORDS (Continue or reverse side if neceseary and identify by block mumber) 


guided projectiles, precision guided munitions, guns, qun Jaunched ramjet, 
soldi fuel ramjets, antiship missile defense, 


| 


20 RUS{EACT (Continun on “aroree side Ht aecoesary end identity by block numbar) 
V/This report covers work done on performance analysis of a 5 inch, 54 
caliber gun-launched guided projectile with solid fuel ramjet (SFRJ). 

A computer program (TRAJET) was developed. The program contains ramjet 
and trajectory analysis. The ramjet part considers conical shock wave losses, 
inlet boundary layer losses, normal shock josses, subsonic diffuser recovery, 
expansion into combustor losses, heat losses at the combustor and nozzle losses 

A flat earth trajectory with drag and thrust was donsidered. The various 


os 


teh BS x Jets 


att 


2 


RR EES 


prey 


a 


OD (jan os 14730 col rion OF | NOV 45513 OBSOLETE UNCLASSIFIED 


SN 112-514-6401 ne ne arr eear ren a TE ae 
SECURITY CLASSIFICATION OF THIS PAGE (Phen Data Entered) 


aoe See ti et eee “Soe eee on 


UNCLASSIFIED 


7 . -- heb Te ee Pan State be ceced- 
20. Abstract continued: 
drag coefficients which were considered are: cowl drag coetficient, skin drag 
- coefficient, wing (cr fin) wave drag coefficient and wing (or fin) friction 
drag coefficient. Base drag is assumed to be zero duc tu the jet from ramjet 
nozzle.\ : 
It was found that the cowl drag coefficient has a major influence on the | 
results,’ Therefore, a separate program (AERO) was developed to calculates this 
parameter. F 
~The 5"/54 solid fuel ramjet has a capability to produce fuel spect fic | 
impulse in the order of 400 - 900 sec. depending mos tly on the flight eit: tude. 
The thrust coefficient varies in the range of 0.3 * 0.1 depending on the 
internal areas. Noe te 
A range in the order of 50 miles can be achieved with the ramjet 
operation compared to only 13 miles achieved by tha conventionas projectile. 
At low-altitude launch, a range of over 18 miles can be reached in Air-Nerense 
Scenario. * The ramjet propelled projectile reaches the ranges mentianed above 
at high Mach numbers (Mo > 1. 8). MIt is, therefore, clear that the ramjet 
concept provides significant imp vement and has an Anti-Ship Missile 
Defense (ASMD) capability. 
| 
| 
\ 


mae 


a 
Oe an 
i H 
ae 
t i : 
i ae 
| 
{ 
______ UNCLASSIFIED __ | | 
ait FET LGTY ashe Stic OF THIS PAGE When Date Freres: t 4 
: | 
[1 wich Semeneterrpmmenr arom mts oo a = ~~+ antennae annneh ith tints: atien tuademesenee * ' j 


= ¢ 


Den Jacohs of the Naval Postgraduate School Educational Media Uepartment. 


Aig ne : a t 


fas RE 2 ee ee tee er ee eR Teter oe 


ACKNOWLEDGMENT. 


The author wishes to express his sincere appreciation to Distinguished 
Professor Allen E. Fuhs for his invaluable guidance, help and friendship. 
The ai'thor would also like to express his gratitude for Frofessor 
T. H. Gawain for his assistance and cooperation in calculating the cowl 
drag coefficient. 


he author also gives his thanks to Professor D.W. Netzer for his 


technical assistance. 

The author would like to thank Mr. Uri Katz for his help and for the 
fruitful discussions we had during his stay as a research associate at the 
Naval Postgraduate School. 


The illustrations and graphic arts for this report were accomplished by 


His help is appreciated. 


. cu brid ifeteit Tim a ee Oe ee ne arte 


epee S 


ee eaion ror 
oan ec anesthe eee > 


ae. 
x =: ne nee = ee. 
Pima eee - Sr ae ae ae mm tte eA lO er a ne ae 


2 


| 
| 


| 


ABSTRACT 

This report covers work done on performance analysis of a 5 inch, 

54 caliber gun-launched guided projectile with solid fuel ramjet (SFRJ). 

A computer program (TRAJET) was developed. The program contains 
ramjet and trajectory analysis. The ramjet part considers conical 
shock wave iosses, inlet boundary layer losses, normal shock losses, 
subsonic diffuser ~ecovery, expansion into combustor lasses, heat losses 
at the combustor and nozzle losses, 

A flat earth trajectory with drag and thrust was considered. The 
various drag coefficients which were considered are: cow! drag coefficient, 
skin drag coefficient, wing (or fin) wave drag coefficient and wing (or 
fin) friction drag coefficient. Base drag is assumed to be zero due 
to the jet from ramjet nozzle. 

It was found that the cowl drag coefficient has a major influence 
on the results. Therefore, a separate program (AERO) was develcped to 
calculate this parameter. 

The 5"/54 solid fuel ramjet has a capability to produce fuel 
specific impulse in the order of 400-900 sec. depending mostly on the 
flight altitude. The thrust coefficient varies in the range of 0.3+0.] 
depending on the internal areas. 

A range in the order of 50 miles can be achieved with the ramjet 
operation compared to oniy 13 miles achieved by the conventional pro- 
jectile. At low-aititude launch, a rang. of over 18 miles can be reached 
with the ramjet version. Launches at high elevation angles can be 
useful in air-defense scenario. Tne ramjet propelled projectile reaches 


the ranges mentioned above at high Mach numbers (M)>1.8). It is, therefore, 
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clear that the ramjet concept provides significant improvement. and has 


an Anti Ship Missile Defense (ASMD) capability. 
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1, INTRODUCTION 
1.1 Background 

This report covers work done on propulsion and flight me- 
chanics of a gun-launched guided projectile. Gun-launched guided pro- 
jectiles have been developed by Martia Marietta for the U.S. Army 
(Copperhead) [1,2] and more recently, also, for the U.S. Navy [3,4] 
(5"/54 Mark 46). The USN round nas solid rocket propulsion. 

The addition of a \iquid fuel ramjet (LFRJ) to the Navy's 
version, was examined in the past by Brown [5]. This report concentrates 
on the addition of a solid fuel ramjet (SFRJ), instead. 

Tt is believed that solid fuel ramjet has some potential 


advantages compared with the liquid fuel. Some of these advantages 


are: 

-Simpie desion 

-High reliability in operation 

-Low cost 

-Fuel control system not needed 
On the other hand, there are alsu a few disadvantages to SFRUJ compared 
to LFRJ. Some of these are: 

~Difficult to control magnitude of chrust 

-Difficulties in achieving high combustion efficiencies 

In both cases, the addition of propulsion improves dramatically, 

the performance of the projectile by multiplication of range and enhanced 
maneuverability. Even more; to operate and nroduce thrust, the ramjet 
engine depends anly on its forward motion at supersonic speeds and does 
not employ any moving parts. This fact, which is especially emphasized 


in SFRJ, Teads to some advantages of the ramjet concept over the other 


> 


propulsion alternatives, at supersonic soeed. On the other hand, the 


supersonic operating regime. The boost required causes some system 


| ramjet engine requires an auxiliary booster to accelerate it te its 
difficulties, But, whife solving these problems, the ramjet system 
| becomes even more attractive for use with gun-launched guided pro- 
jectile, like the U.S. Navy 5"/54, 

A computer program was developed to analyze the performance i : 
| of the SFRJ, The computer program was writtan for the IBM-370 
computer at the Naval Postgraduate School, Monterey, California. HTPB | 
| was selected as a fuel, but performance with any other fuel can be tested, ' 
; using the same model. A flat earth trajectory with drag and thrust was | 
1 considered. Using solid fuel, a thrust-equal-drag trajectory is more | 


‘ difficult to achieve with SFRJ. Therefore, most of the results given 


in this renort eliminate this case, and the exact value for drag, at ; - 
Poy 


each point, was calculated. However, if cesirad, it is possibile to 


: fn 
change the air mass flow, in order to obtain thrust-equal-drag flight. tee 
The computer program can calculate this case also. | 


1.2 Basic Concept 


! The ramjet engine consists of an air inlet, which serves as j 
4 - a diffuser, a combustion chamber, and an exhaust nozzle [6,71. The | 
diffuser admits air to the engine, which is mixed with fuel (solid or 
liquid) at the combistor. After the burning process, which adds heat 
to the flowing air within the system, the gases are transfered to the 
j 
| 


nozzle. The nozzle ccnverts part of the thermal energy into kinetic 


The areas inside the ramjet engine are usually divided into 


be energy to produce thrust. 
| 
{ 
s 


six Stations, as illustrated in Figure 1.1. Station 0 defines the cross- 
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section area of the stream tube captured by the inlet of the projectile. 
Stations 1 and 2 identify the diffuser. Station 1 itself is at the 
throat of the diffuser, but sub-station 1C is after the conical wave; 
sub-station 11 and Sl are ahead of the normal shock wave, located at 
the throat, or at the actual place, respectively; sub-stations 12 and 
$2 are as above, but behind the normal shock wave. Stations 3 and 4 
refer tq the entrance and to the exit of the fuel grain within the 
combustor, respectively. Note that both A3 and Ay increase with time 
as fuel burns. Station 5 and 6 belong to the nozzle's throat and to 
the exit of the nozzle, respectively. 
1.3 Data 
1.3.1 Dimensions 
In order to be compatible with the Navy's 5 inch, 54 caliber, 
Mark 46 gun mount, aS modified for gun launched guided projectiles 
(Figure 1.2), a set of requirements were adapted initialty: These were: 
a. External shape of existing 5” guided SAL projectile. 
b. Length ~ 60,97" 
c. Length of combustion chamber - 23" 
d. Total weight - 104.7 1b. 
e. Muzzle velocity - 2500 ft/sec. 
Typical values for the internal areas in the ramjet within the 
projectile are (units - 3q. in.): 
A. Aq Ay Ay Ag As 
13.3 5.2 2.6 4.3 8.2 7.5 


Refer to Figure 1.1 for definition of the areas. A, is a reference area 


For typical flight Mach number of: Mo = 3, the appropriate Mach numbers 


at the main stations are typically: 
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1.3.2 Combustion Process 
Some of the losses in the total pressure were taken to be 
constant. These are: 
a. Inlet boundary layer losses (t'y3 typically = 0.93). 
b. Subsonic diffuser recovery (n'y: typically = 0.93). 
c. Nozzle losses (m3 typically = 0.96). 
These values are typical and are not expected to vary too much. 
All the other losses in total pressure, which are dominant to 
the projectile performance,tvere calculated. These are: 
a. Conical shock wave losses (tr,)- 


b. Normal Shock losses (Ts). 


c. Losses due to expansion into the combustion chamber (n,). 


d. Heat losses in the combustion chamber (n,). 
See section A 1.4.2 for definition of various 7 . Combustion efficiency 
was taken constant: ny = 0,90 
Air heat capacity ratio was also taken constant: tg" 1.4; however, the 
value for the gas heat capacities ratio of the combustion products (vp) 
‘was calculated from thermodynamic data for Hydroxy Terminated Polybutadiene, 
HTPB, burned in air [8]. The stoichiometric chemical reaction of HTPB 
burning with oxygen is as follows: 


4 Cag) 93% + 393 OQ + 292 CO, + 206 H50 (1.3.1) 


1.3.3 Trajectory 
In the trajectory part of the program, the various drag coef- 


ficients were calculated. Those are: 


a. Cowl drag coefficient (Con) » [9-12]. 


e 


4 


a 


Newt 


| 
i 
1 


f 
t 


r ~ — Sn 


an 


~~ 


b. Skin drag coefficient (Cys) [13-15]. 
c. Wing (or fin) wave drag coefficient (Corny) > [16 - 20]. 

d. Wing (or fin) friction drag coefficient (Coup) > (13, 19]. 
Base drag (13, 21-23] is assumed to be negligible due to the jet fr n ramjet 
nozzle. 

The model, which was chosen to calculate the cowl (nose) drag 
coefficient, was: based on a theoretical development done by T. H. Gawain 
(9]. The modified program (AERO) is listed in Appendix G. However, 
program AERO as it is, appears to be too long to be used directly in the 
iain program (TRAJET). Therefore, best fit curves for calculated results 
from AERO were used in TRAJET. 

Skin drag coefficients were calculated for either laminar or 
turbutont flows. The same routine was also used to caiculate wing or fin 
friction drag coefficients. To calculate the wing wave drag coefficients, 
a psuedo 3-dimensional model was developed. 

The program alsc has an option to calculate a trajectory of a 
projectile without propulsion. In this case, the drag coefficients which 
are calculated are: 

a. Nose drag coefficient (a different model than the above). 
b. Base drag coefficient 
c. Skin drag coefficient (as above). 

1.3.4 Air Nefense Scenario 

In the air defense scenario, the program takes into account only 
cases in which the projectile exceeds a Mach number of at least 1.8. This 
value of minimum Mach number (XM,) can easily be changed. 

1.4 Results — : 
a Each section of the program was developed, tested and run sep- 


arately. The ramjet part was first run without the trajectory part using 
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constant altitude (Typically - 10,000 ft). The same was done with the 


trajectory part using vacuum case, thrust-equal-drag flight, or constant 


thrust case. The final version was programmed to give the fcllowing 


optional printings: 
a. Loop on all possible values af Ad/A, and Ax/A., and print 


summary tables only. 
Print detailed, time dependent, tables for any specific area 


ratio chosen: 
- Results from combustion process (file name: CMB D) 


- Results from trajectory process (file name: TRJ D) 


~ Various drag coefficients (file name: ORG D) 


c. Detailed print of every step during the calculation, for 


checkup. 


d. Variation of the above: 
Detailed print of cases that were found not to be 


suitable: 
- Reasons only 


- Full detailed parameters 


on Mach numbers, also (output of subroutine 


' 
- Loop 
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2. TYPICAL RESULTS 


Figure 2.1 presents the dependence of the fuel specific impulse (Teo) 


in sec.) of the ramjet on the projectile Mach number at various altitudes. 


~ = ae SS 


It appears that the 5"/54 ramjet has a capability to produce fuel specific 
impulse in the order of 400 ~ 900 sec., depending mostly on the flight al- 
titude. The dependence of Tey on the flight Mach number is weak. 
Figures 2.2 and 2.3 present the dependence of the thrust coefficient 
(Ce) en the internal area ratios Ad/A,. and Ac/A... In figure 2.2, the change 


of Ce with altitude and with Mach number is also presented. The thrust coef- 


ficient (Ce) varies in the range of 0,3+0.1 while Ad/AL. changes from 0.25 to 
0.40 and As/A,, changes from 0.42 to 0.26, 

: The correlation between the fuel specific impulse (T,)) and the thrust 

coefficient (Ce) is presented in figures 2.4 and 2.5. in both figures, a Mach ae 

number of Mo = 3.0 was selected. In figure 2.4, the correlation was checked 


at various altitudes and at various Ad/A; area ratios. In figure 2.5, various 


internal area ratios (Ap/A,» A./A..) are presented. 
= or More detailed results are presented in Appendix H. The dependence of the 


! projectile performance on the other internal area ratios (A, /Ag; Ao/Ags A3/A,) 


1! ' _ Was also checked. Some typical results are presented in that Appendix. 
An altitude ys range dependence for various elevation angles is presented 


in figure 2.5. A range of over 80 km can be achieved with the ramjet operation, 


i jectile. A low-altitude launch (in this figure, an elevation angle of 15° 
was selected) is also presented reaching a range of over 30 km with the ramjet 


| 

compared to only slightly mere than 20 km achieved by the conventional pro- | 
| 
operation. The high elevation angles are mostly used in air-defense scenario. | 
| 


The drag of the projectile when the ramjet is not operating, for example, 
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after burnout, was not determined. The computer program does not account 
for the drag increase due to the ramjet not operating. Consequently some 
of the trajectory curves in figure 2.6 are in error, However, for tra- 
jJectories at low gun elevation, the ramjet burns all the way to splash. 
These trajectories are accurate. For all trajectories, the curves are 
accurate to the peint of ramjet burnout. The trajectories of interest to 
air defense are accurately calculated. Trajectory of thrust-equal-drag 
(vacuum) case is shown for comparison. 

Results for air defense scenario are presented in figure 2.7. Only 
ranges where the projectile Mach number exceeds at least 1.8 where con- 
sidered. The area ratios where chosen as specified in the figure. The 
two cowl angles, shown in figure 3.1 were 20° and 9.5° respectively. 

The gun elevation angle was varied from 7° to 80°. The change of atmospheric 
conditions with the altitude was taken into account. In table 2.1, some 


typical results for "Surface-to-Surface Mission" are presented. 
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TABLE 2.1 
Solid Fuel Ramjet Propelled 5"/54 Projectile 
Surface-to-Surface Mission 
Ranges (km) vs Gun Elevation Angles 

Elevation 
Angle ; 7 25° 45° 65° 80° 
a. Ramjet?) 80.3 15.9 5.6 
b. Projectile 9.2 7.5 20.2 16.3 7.6 


fabinee 


Without 


Propulsion 


Note: 1. Area ratios as in Figure 2.7 
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3. Discussions 

Looking back at figure 2.1, the dependence of Tey on altitude and on 
| Mach number should be explained. We shall do that by using the equations 
described in Appendix A. 
From equations: (1.4.9a), (1.4.12), (1.4.16) 
together with equations: (1.1.1), (1.1.2), (1.1.3) 


one obtains: 


L ley = Sp XX (3.1) 
| where: gh Wet (3.2) 

y and: Cp = ky > ky (3.3) 
be where: Xo = ky “1X11 ~ (3.4) 
X1 = kKgMglitkc(Py fig) | (3.5) 


The parameters ky to Ks are functions of the various area ratios, the 


Me Nee Reamer meee tt 


temperature of air (T)), the heat capacity ratio of air (ya) and the 
perfect gas constant (R,). These parameters are assumed to be constants 


hot, in discussing the influence of the change in altitude and in Mach number 


| on the vaiue of Teo The altitude dependence is mainly due to change in 


0.4 


dependence of Iep on pressure is approximately Po 


| 

\ 

eos 
fe ~ atmospheric pressure (Po). In the conditions chosen for figure 2.1, the 

at 45 =3. That means | 

{ 


that, in the region mentioned, Ion pressure dependence ts mostly due to | 


change of X, (equations 3.1 & 3.2). From the same equations, the Mach 


Li number dependence of Teo can also be explaines. At high altitude, the 
. change of Ks (Wr) is very close toa the Ce dependence on Mach number 
|: (ig) and therefore 1.5 is almost constant while changing My. On the | 


7h 
other hand, at sea level, the change in a(t 28 is smaller than that of 


i 2 
Ce(My 2V4y and therefore Ten changes with My as shown in figure 2.1. i | 
| 19 
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Figures 2.2 - 2.5 present similar dependences of the ramjet performance, 
and can well be understood using the same equations. 

Testing these results, together with those presented in Appendix H, 
the design of the ramjet internal area ratios can be completed. The 


results are listed in Table 3.1. 
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Table 3.1: Ramjet Design ! 

Dimensions 

Externai diameter = 5" ; 

: " Total length = 60,97" 

| Total weight = 104.7 Ib (47.5 kg) Ps 

! Area Ratios: | 

7 Ag/A, —Ay/Bg g/g a/A, g/g AY | 

. 0.25 0.47 0.887 0.426 —0.26 1 - 

j Reference Area: i 

a A, 219.3 sq. in. (124.5 on”) | 

: : Combustor ao 

| : Solid fuel: Hydroxy Terminated Polybutadiene (HTPB). ly 

Fuel weight: 3 kg | | ; 

qe Fuel density: 971.56 kg/m? i 

: | Fuel specific impulse (1): 400 - 900 sec. | 

val Booster an | 

iy Booster weight: 2 kg I 

: “Booster density: 1650 kg/m? | ; 
Z Booster specific impulse (1.5): 240 sec [+ i 
: : Performance | 
} Muzzle velocity: 762 m/sec | . 
* "Velocity after booster: 863 m/sec | 
DS Thrust Coefficient (Cp): 0.2 - 0.4 | 

1. 

; ' 
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In figures 3.1 - 3.4, the designed ramjet concept is presented. The 
guidance and the control sections as well as the warhead were not redesigned, 
The location of the tailfin is described in figures 3.3 - 3.4 [White,24]. 

The configuration of the solid fuel ramjet presented here, is in accordance 
with the design of the liquid fuel ramjet done previously by Brown [5]. 
The new design is also in agreement with the Navy's requirement to be 


compatible with its 5"/54 Mark 46 gun mount, as modified for gun launched 


| 

guided pr:jectiles. 

be 4 Figure 2.6 presents a full-range comparison of the SFRJ 5"/54 pro- 

| jectile performance with that of the conventional projectile. The 

improvement in performance of the ramjet-propelled, cuided projectile is 
very significant. The extended range of the ramjet concept can be used in 


a "Surface-to-Surtace Mission® (Tabie 2.7). The main improvement of the 


Ace 


! 

| | ramjet concept might be in an "Air-Defense Mission" (figure 2.7). A 
minimum Mach number of Mo = 1,8 was assumed for both the ramjet concept 
and the conventional projectile. At lower altitudes, the projectile is 

: limited by burn-out of the fuel. At higher altitudes, the performance 
boundary is due to ramjet unstart. The corventional projectile is limited 
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| 
| ‘ roe Sale, Su P ~ sem de * 4 i 
a by the Mach nuiiber cecay to values Tess than 1.8. It ts self-evident 
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ore eee 


that tne ramjet concept provides significant improvenient, and, therefore, 
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has significant ASMD Capability. 
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Figure 3.2 Solid Fuel Ramjet Propetied 5"/54 : | 
Projectile: Design, Section B-B 
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4. Conclusion 

Based on the computer model discussed above, the ramjet-propelled, 
guided projectile provides significant improvement over the conventional 
projectile. The fue’ specific impulse of the SFRJ is in the order of 600 
+200 sec, depending mainly on the altitude of the projectile. The 
appropriate value of the rocket is only 300 sec. The thrust coefficient 
varies from 0.2 to 0.4, depending on the atmospheric conditions (altitude) 
and on the geometry of the projectile (internal areas). Therefore, the 
ramjet-propelled, guided projectile reaches a range of about 80 km compared 
to range of slightly more than 20 km in the conventional projectile. The 
improvement of the ramjet concept might be in both “Surface-to-Surface 
Mission’ and in “Air-Oefense Mission”. It provides an ASMO weapon which 
is compiencntary to guided missiles. 

More work is required toe design an optical system incorporated into 
the inlet. For satisfactory ramjet performance at the flight Mach numbers, 


the lens must be a conical shape. 
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Appendix A: SOLID FUEL RAMJET: EQUATIONS 
Al. Combustion 
Al.1 Computation of fuel - regression rate, weight rate of 


burning fuel, and fuel - air ratio 
Define G as the weight flow rate of air per unit through the 


| entrance port to the combustor; see figure 1.1. Uence G is given by: 


a 
G= G (1.1.1) 
and has dimensions of Tb/sec.in’. The simple form of the regression 


rate of a solid fuel, *, iS given by: 
> = ag" (1.1.2) 


Where a and n are empirically determined constants. The dimensicns for 


Se 


Li f are in/sec. Knowledge of the temperature dependence of the regression 
ie rate will allow the use of a more accurate model instead of equation 


; (151.2); 
i The weiqht of fuel burned per unit of time is as follows: 


We = Pe MrOeb3 (T4d2) 


we a a MNO EE 


BE Where p~ is the density of the fuel in Ib/in’, D, is the inside dia- 


rel 


. meter of the fuel grain. Note that 0, increases as the fuel burns. 


. The length of fuel grain is Ls inches. For the case of HIPB, the value 


of 0.035] ib/in? was taken. In an actual solid fuel ramjet, * varies 


4 


: alona the grain; * is largest in the region immediately downstream of 


a 


the entrance port of the combustor. However, for the ramjet model 


+! \ : 
ne ! developed here, the value cf * is assumed to be constant along the grain. 


iahsdatethashanetlionton sean tahe erienetedan at ae 


ra : Consequently: 


t 
D3 = Dag + 2) : tdt (1.17.4) 


The integral form for the change in grain internal diame - is used 
since * may vary with time. The initial grain inside diameter is D2 
inches, and the inside area for fuel grain as a function of time can be 


: ; written as: 
| [fF agp gees 
a Ag = Z . 7 439 + Se rat] (1.1.5) 


_ By definition, the fuel - air ratio is: 
Ps Pet (1.1.6) 
The value for We is obtained from weight flow through the inlet, and the 
value for We is calculated using equation (1.1.3). The total mass flow 
ey through the nozzle is given by: 
dee Wy We +R, (1 + FD City 
For HTPB burning in air, the stoichiometric value for f js 0.0728; high 
i combustion efficiency is difficult when f is less than 0.025. 
! Al.2 Computation of Combustion Exit Condition 
| Combustor exit conditions are specified by four quantities as 
| follows: stagnation temperature, Trq> Ops stagnation pressure, Porge 
psi; ratio of heat capacities, Yes and gas constant, Re in/ 9p. In 
the computer program, the appropriate mks units are used, i.a. Oy. kgins 
my, respectively. To determine the exit conditions, two input quan- 
bet tities are needed. These are fuel - air ratio and stagnation temperature 
at the combustor inlet, Tyo: Note that Tr3 = Trp has been zssumed. 

From the thermodynamic data for HTPB burning in air, one determines 
Tyg (th), Yer Re. The symbol Tyq(th) is a theoretical temperature which 
results from 100% combustion efficiency. Introducing the definition of 


combustion efficiency yields: 


ee en ee en See Soe 


Tie 
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As discussed previously, a constant value of Ny equal to 0.9 has been 


assumed. 
The value for Pra is calculated based on one-dimensional, choked 
nozzle flow. A certain value of Pra is required to force a certain weight 


flow, Wy, through the nozzle. Assume that V¢ remains fixed through the 


nozzle. Define a function of Y¢ as: 


3 \Y¢ + 1)/{2(v¢ -1)] 
pa Fees | (1.2.2) 
Five +1 
Define a characteristic nozzle velocity, 
c*, m/sec: 
VG Re tT 
c* = 7 — (1,2.3) 


where g is the acceleration of gravity and has value of 9.807 m/sec’, 


The required value for Pra is givan by: o 


Wyc* 


Prq * Ke (1.2.4) 


The decrease of flight stagnation pressure, Pro» by inlet and combustor 


losses must not be too large. If the inlet and combustor do not provide 


the required Pra» the inlet will unstart and Ws will decrease. 
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Al.3 Computation of Nozzle Exit Conditions 


The relation between the area ratios and the Mach number is well 
known by the formula: 
| ree (ve + 1)/[2y¢ - 1)] 
i = Me (1.3.1) 
+ tee Saige 
6 


Ags Ae are the areas at the throat and at the exit of the nozzle, re- 
spectively. Knowing y-, the exit Mach number (Me) can be calculated 


. for any nozzle area ratio (Ay/Ag). This indirect calculation is done 


‘ in subroutine CALCM, using Newton - Raphson's iteration routine. 
t- 
tT: The iutal pressure at the exit of the nozzle (py_) is defined by: 
“ao, 
: Pre = Pry, (1.3.2) 
ae where the total pressure at the exit of the combustor (p_,) was calculated 
i! 
i previously (e.g. 1.2.4). 
: Knowing the total pressure at the exit of the nozzle (Pr¢) and the 
| 
ou Mach number at this point (Me). the exit pressure (D¢) can be calculated: 
- ie Z < “a 
ee : Pe Pre 1 tan MG es 
ee 
be Al.4 Computation of Thrust and Thrust Coefficient 
a oo 
a6 A).4.1 Thrust Coefficient (Ce) 
qo! OS 
The thrust of the engine is the net rate of change in 
momentum at a steady state condition, and is given by: 


Fos Pehg + hep - Porg - Uy - Po (Ag - Ay) + Pa lAg-Ag) 
(1.4.17) 


= Peds + MgUe ~ Podg - MU 


Ga aera 


i: SIs ce ere _— — Rs oe 
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where Uys Ue» tho: the are the velocities and mass flow at the inlet 


entrance and at the nozzle exit, respectively. 


From the continuity equation, the foliowing relation arrives: 


mU = pA 


(1.4.2) 


Substituting for the density (o) from the perfect gas equation of state: 


ais Rr (1.4.3) 
gives: mu = By ura (1.4.4) 
From the definition of Mach number and speed of sound: 
ue = wae ; a? = yRT (1.4.5) 
Therefore: nh = Ben? yRTA 
nl = oM@yA (1.4.7) 
co Substituting (1.4.7) into (1.4.1) gives: 
a A. A Jf 
i F = peA(l + Y¢te?) PoAo(x. + yaMy°) (1.4.8) 
The thrust coefficient is defined : 
7 pounce 
. f "af, (1.4.9) 
et where: dy = A0qUy” = bY4PoMye (1.4.10) 
: i 0 0 a 0 
{ 
Combining equation (1.4.9) and 1.4.10) gives: 
i F , 
C, = aera rs (1.4.11) 
asl ee Ply A 
‘ [ x rT 


Substituting for the thrust from equation (1.4.8) turns equation (1.4.11) 


into: 

I , 2Ag/A,. Races 9 | 2Ao 

i. €, 2 (1 + ¥4t.S)-1 | 7 (1.4.12) 
4 Fo yatty? -LPre”Ps #' i 
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Al,4.2 Prassure Losses 


The pressure ratios in the above formula: | 
Pre/Pro 
Py ¢/Pg ! 
can be substituted by a function of pressure losses across the ramjet: ae 
Pre _ Pre Pra Pr3 Pre Pro (1.4.13) | 
18 - 110 2th oly ltée liv 
Po Pra Pr3 Pro Pro Po | 


We define the pressure losses as follows: 


| 
| 
Prg s s N ] 1 s 
Brg T, = Nozzle losses | 
is 
Pra . ; i 
—* = 7 = Ravleiah flow losses ee 
Prg a | oe 
Pr3 _ , Bo 
Bro Te = Combustor expansion losses tag 
T2 ae 
% Hi gig 
Pre me 
me)” Inlet losses = (conical wave loss )*(boundary tayer loss)* ee 
TO ! i 
*( normal shock loss }*(subsonic diffuser recovery) = TET 'TyeT py" Ly | 
Therefore: 
2H 3 
(ng) oo 
p p I(n. Pro pee 
T6 ; » *TO _ 1 a 
Do - Th Th Te uy Th TNS Th Po = be = losses Pg (1.4.14) | | 
Finally, substituting for total pressure ratios (pr,/p,; Pre/p,) gives: 
TO’ FO? "T6/ "6 ane 
Val (ygr1) ! 
fb + Sara? | 
+ 
Pr6/Po _[ittn.) Ol yl) 15) . 
ee j Yer oY ¢/ (Ye (1.4. } 
Prg/Pg  Li=tosses [; sat Z| 
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This relation can be used in equation (1.4.12) which calculates the 


thrust coefficient of the system. 


A1.4.3 Computation of Thrust 


The thrust can easily be calculated from the thrust coefficient, 


using eq. (1.4.9): 

Foz Ce I A. (1.4.9a) 
The dimensions of F are Newtong (after multiplying eq. (1.4.9a) by 
The fuel specific impulse, Top? in 


the acceleration of aravity, 9). 


N/kg/sec, is defined by: 


= F/We (7.4.16) 


Tsy 
is given by: 


The specific fuel consumption, SFC, in kg/hour/N 


SFC = 3600/1, (1.4.17) 
Ramjet performance is specified in tims of the perfcermance paremeters 


and SFC. 


Ce, F, Tey 
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AZ. Check for Choked Nozzle 
The total pressure at the throat of the nozzle is given by: 
Prs = Pra J" (2.1) 


Again, Prg> is the total pressure at the exit of the combustor, and is 


calculated by eq. (1.2.4). un is the nozzle loss. 


It follows that: 
. 2 \rel (yet) (2.2) 
Pp = Pts (Ye 


f 


The pressure at the throat of the nozzle (Px) should be equal or greater 


than the atmosperic pressure (Py). 


Ps> Po (2.3) 
When inequality (2.3) is satisfied, the nozzle is choked. 
A3. Heat Losses at the Combustor 
A3.1 Mach Number 
a. Continuity 

At the combustion chamber, fuel is added and the continuity 
equation is: . 

03 Uz Ag (1 + f/a) = % Uy Ay (3.1.1) 
where Ay and A, refer to the entrance and to the exit of the combustor, 
respectively. By assuming that A, = Ay. the continuity equation 


{eq. 3.1.7) can be written as follows: 


03 Uy 
— (1+ f/a) = (3.1.2) 
4 U3 


Replacing the velocities (U5 and Ua) by the appropriate Mach numbers 
(eq. 1.4.5), turns eq. (3.1.3) into: 


p MA\ iT, ve R 
3 4); 4 '¢ “f 
-2(1 + f/a) { ) oe : (3.1.3) 
4 MSA) 73 Ya % 
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| 
Where Rae Re are the gas constants of air and of the combustion products, 
] 


respectively. Hence: 


: 6 
(*2) (2+) ( rt) Ee 1 
| Mgh Npg/ \T37 WAR 1 + 7a (3.1.4) 
b. Momentum 
; Applying the conservation law of momentum to the discussed 
problem: ; 
| ; 
i Ps + P4 us" = Py + Pa uv,” (3.1.5) 
t 
| From the definition of Mach number (eq. 1.4.5) and the perfect gas 
equation of state (eq, 1.4.3): 
te 2 xs 2 
4 0, U" = ¥; py My (3.1.6) 
jj Substituting equation (3.1.6) in equation (3.1.5): 
4 2 
| Ph ge a (3.1.7) 
4 Va ye WS 
i 4 Y¢ 4 ! 
i ! i 
Again, from the perfect gas equation of state: i 
J 
a | Pq [Pa (4) (3.1.8) 
a 3 (03 /\T3 
Ya 
aie Substitution of (3.1.8) into (3.1.7): 
eee 
| | 
G) : (2) 1+ y,M, (3.1.9) 
. 2 a tse ! 
i 
ve | 
i] Substituting (eq. 3.1.9) turns equation (3.1.4) into: | 
z 35 2 i : 
ie (e) (*) . flit vals \ [tee 1 
‘| LANEY ae vey Yara T¥ 7a (3.1.10) 
ie 
| 
ij 
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Replacing the temperatures T3 and ‘4 by the appropriate total temperatures: 


<1 
Tas Tes Yj 2) 
i ni + ten, 
and assuming that the change in total temperature up to the entrance of 


the combustor, is negligible (Tr = Tr3)> one obtains: 
a 


45 2 Yel 2 
«{' + anes 3 \ pit See Yee 1 447) 
2) (et y_R 1 + f/a 


1+ Myo i. tay,2 aa 
ies 
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c. Solution 
The Mach number at the exit of the combustor (My) can be 


solved knowing the conditions at the throat of the nozzle: 


Y¢r 2 (3.1.12) 


The computation is again indirect, using Subroutine CALCM. Knowing 
My» Equation (3.1.11) is usad to compute Ma. The solution ts received 
LS da fenn the 


by iteration. As first approximation, hig can be soived from 


following equation: 


( ras} = (3.1.1a) 
M 
; : ; 2 P + 
where: ae 1} + eM, (=) Yara qq ve f/a) (3.1.13) 
4 Tyo vere 
consequentTy Ce aoe 
+B- Bo -4 
Rees 
May = ast (3.1.14) 


“The subscript N shows that the computation was made from the nozzle 


37 


1) eee lt RETIN SIPEG RCN: REE ROE oo ene a ee 


TES tne a le al 
pe 
7 a seme a a pr re aa SaaS Peers si SE is tae 
woae “> Aorta MRT Oy emer ee toa 


he a enema, soot . 
tae = ED are ne epreen . es 
— ec - : 
i eee ee Soiree eee Set 
Pie Aaaaens 


| direction. Now, B is changed to be: 


! 
pea 3 
a 2 
1+ M. } 
oe eae —_— (3.1.15) 
(4 ny? 
Substituting 


j 

t This expression is consistent with equation (3.1.11). 
{ 

back into (eq. 3.1.14) gives an improved value for M3. This procedure 


can be repeated several times, but it was found that even after two 


interations, the value received for Moy is accurate enough, due to the 


small change in B resulting from (eq. 3.1.15). 
A3.2, Total Pressure 

de d By definition: 

{ 

“ie vel lyer}) 

‘ Pr 7 yz ‘ } FONE 

ie is ee (3.2.1) 
“| i ‘ 

j Consequently: . 

4 a al (rat?) 

| Prs = 53 Poe bd “s (3.2.2) 
me Pra Py i" Ys* 

Rg 
at a 

R ‘ Substituting for (P4/P4) from equation (3.1.7) resuits: 

i ) ; Wai Yq/ (vat!) 

: Pran ~ Pra coed ve a 

ie a) ee ene? (3.2.3) 


Where the subscript N was defined previously. 
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A4. Computation of Mach Number and of Total Pressure at the Various 
Stations of the Inlet 


Ad. Initial Conditions 
In the previous sections (Al - A3) the pressure conditions at 


the combustor and at the nozzie region were calculated. Here, the 


pressure conditions at the inlet will be calculated independently. 


Knowing the total pressure conditions at the various stations of the 


{ 

| 

| 

: inlet will allow the creck of whether the inlet can supply the amount 


l 
of air needed by the combustor. As will be seen afterwards, this check 


will also allow to specify the location of the normal shock wave at the 
a inlet. 

Aoi 

| Assuming that the static pressure (no). the static temperature 
(Tp) and the flight Mach number (Hq) are known from the trajectory part 


of the program, the total pressure and the total temperature can be 


oe eee ee 


But for flight Mach 


| 
| 
calculated: 
| val bg) 
ts v.> Y 
joy a Ya a” 
. Prg = Pg (1 + 3,73 (4.1.1) 
ad T ra 42) (4.1.2) t 
aS trot to + ne : [. 
| The weight flow through the inlet is given by: ' 
a W. = 6p Un Ay pe ’ 
4 a * 9 Yo Ao (45) i? 
Fel Ag 
j i.) Usually, when flight Mach number is equal or greater than the inlet “ 


design Mach number, the value for Ac/Ay is unity. 
A 


et 


The additive 


te ee, 


{ 

{ 

i 

bo ae : . 

(G3 number less than design Mach number, Ac/Ay becomes less thar 1.0. 
i i 

ie value of 0.9 was selected as a constant value for A;/A,,. 


drag due to Ac/Ag < | was ignored. 
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* area behind the conical wave (Pro. Myce Arc» respectively). 
A4.2.1 Pressure 
The pressure coefficient can be defined as follows: 
Dy n7P 
Co ee OO as (4.2.1) 
p Va p M 
> 8) 0 
For a cone, the pressure coefficient (C,) can approximately be for- 
mulated as: 
- 1.69 
: 0.096 | fa 
Cy | 0-083 + ae ] (#5) (4.2.2) 
0 
where a is the cone half angle. The difference in pressure on surface 
and behind shock wave is ignored in this model. 
Knowi na Cy the pressure ratio (Py ¢/Pq) can be calculated (4.2.1) 
D Y 
1c ay 2 
— =1t M 2. 
The Same pressure ratio, is alsa related to the Mach number, normal to 
the conical shuck wave (mM). Using (ea. 2,48a) in ref. [16], one can 
get: 
Pie oy a NON) (4.2.3) 
Pg Ya"! 
Knowing (Py ¢/Pq) from equation (4.2.1a), M, can be calculated from 
equation (4.2.3) : 
Pac v4" 3 (4.2.3a) 
Mo o= i+ (¢ 7 1) = 
n Po “Ya 
40 
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A4.2 Conical Shock Wave Loss 
In this section, the conical shock wave loss will be computed; 


the calculation results include the total pressure, Mach number and 


1 Auingcaihd ou ieee vuabet Stee, a 
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After computing the pressure ratio due to conical shock wave (equation 
4,2.la) and the Mach number normal to the cone (equation 4.2.32), one 
can use equation (2.54) in reference [16] to compute the total pressure 


ratio at the conical shock wave: 

- Tee 2 4 y,/(v,-1) 

a x Pric 2 cr zy, 2 | Wy, 1) Yat) M atta 

Ll aie Ei Biers g (M °-1) (4.2.4) 
| Pr) L Ya aaa Yat! M42 

| From equation (4.2.4) one obtains: 

| 


— Pric 7 Pra 7% (4.2.4a) 


t : A4.2.2 Mach Number Downstream of Conical Shock Wave 


Using figures 4.1 and 4.2, the wave angle (8) can be defined 

s! as: M 

0 = arcsin i) (4.2.5) 
u 


On the other hand, the deflection angle (@) is defined from equation 


(4,10) in 1 oference [16]: 


(m_7=1) 
iu ] (4.2.6) 


@ = arctan | zeot (8) a ae eae 
My {y,tcos (28) }+2 


The Mach number behind the conical shock wave (Mic) may, therefore, be 


ee - obtained using equation (4.7) in reference [16]: 
a . y.71 
: 1+$—4* 44 
te = n Pea 
: Me? | = er eee | ee. 
(3 U2 a 2 
me PRE. 


The relation between the deflection angle (4) and the wave angle (8) for 
various Mach numbers are shown in figure 4.3 (reproduced from figure 4.2 


i in veference {16]). 
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A4.2.3 Area Ratio 
j 
Using figure 4.2, the area ratio behind the conical shock 


wave can be obtained: 


IE Cer apes) 


1¢ , Letan(e) case) 14.2.8) | 
{ 


| 

! 

Equation (4.2.8) gives the area normal to the flow at the inlet lip. 

| A4.3 Boundary Layer Loss | 
| | A4.3.1 Mach Number | 
| The area ratio is related to the appropriate Mach number by | 
| 


! by the formula: 
ei : 
A (y,' ae 


4 ; ry 
a Be" ac 2 
e me? 


a 
( 
) 


Where Ay 
yr, FU) /120yg-1)) 


* is the area at the throat of the inlet. Similarly, 
a 
os (4.3.2) 


iy and M., relates to the area and to the Mach number ahead 


Ll aw, A 


a, } AnSTe “TH 
of the normal shock wave. Dividing. these two equations gives: 
(ygti)/f2(y,-3)] 

(4.3.3) 


we ite 
£ 1 
he Jas 


(4.5.2) 
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a | 
| | 
Knowing Mics Ya> Ayes Aqy: equation (4.3.3) can be used to calculate Maye ! 
| indirectly, by subroutine CALM, whicn was mentioned previously. Mid 
should be supersonic (My) > 1) to prevent unstart conditions. ! 
i 
A4.3.2 Pressure : | 
The total pressure in front of the normal sheck wave (Pry) is 
i 
me received from the connection: i | 
os i 
Co Prin Pric (4.3.4) , | 
| Where m)', the boundary layer loss is assumed to be 0.93. 
i 
| A4.4 Normal Shock Loss j | 
i 
The Mach number, behind the normal shock wave is defined as: } ! 
{ 2 ! 
| My? + ¥e-T i | 
Mi2 = 11 a aes (4.4.1) ' 
——~r Ma = 1 
|. ! \Ygr7 il | 
{ The total pressure behind the normal shock wave is defined as: | 
: | 
a4 i 
ie vqth ' ¥a/(¥,e1) L/(y,-1) | 
“| ~ ~ M ) 2y y,-l 
p =p eT a 2 Ya Ee 
| M2 STH perera ar “17 yt (4.4.2) | 
| ohh fer. | . 
| A4.$ Subsonic Ditfuser Recovery | 
“a Similay to equation (4.3.3) ene can odtain: | “ / 
ln ; : ae pein) | ay 
| Bie gg MEN ee | 
; re Ny Y.-1 (4.5.1) 
le 2 a 2 
te 1+ =f M 
a | é 12 Ea 
L le Knowing May? Yo? Ayo: and No» the Mach number at the exit of the inlet | 
14 | 
dj can be comsuted using subreutine CALCS. The total pressure at this station : 
i 
aA is defined as: 


Pro * Prig ™p” 
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Where the subsonic diffuser recovery (1t5") is assumed to be 0.93. 


ye A4.6 Expansion Loss 
: j 
A4.6.1 Mach Number 
} 

| 


On their way to the combustor, the gases coming from the inlet 
expand at station 3; a sudden change in area from Ay to Ay occurs. The 


sudden change in area acts as a flameholder by creating a hot recirculation 


[ 
j 
f 
| 
region. In this section, the loss in total pressure due to this expansion 
is calculated. | 
From the continuity equation: H 

i} 


| PoUoAy + Pyll,A, (1.4.22) 


i From perfect aas equation of state: 


o= Be (1.4.3) | 
| 
| 


ns 


ia. 4 


and from the definition of Mach number and speed of sound: 
M = U/a;: a= VykT (1.4.5) 


Equation (1.4.2a) turns, therefore, into the form: 


; is P3 
Rr, Me Mighty A * ts 3 WyaRT3” Az (4.6.1) 


f 


For sudden expansion of an incompressible fluid, the change in static 


pressure going from small area Ay to large area Ay is given by (16, 25]: 


7 


e + ok 
OP aa es wah a a ee 


By ~ Po 


Where Aas = Ag/Ag. Also, for the incompressible case, the stagnation 


. 
eran % me =. 
1 
eo? 
5 ‘ 


7 re pressura ratio is given by [16, 25]: 
I : 
"toe p G5/P 
ity 7] oe ae 7 2 4 

ai Pr 1 Y¥q5/P (1 Any) { .6.3) 

i> ote 

tj fccordirg 10 equation (4.6.3) as My decreases Pr3 approaches Pro: The 
iz 
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model for sucden expansion with compressible flow is inch more complicated. ‘ | 
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i 
Tt was assumed that static pressure is constant in the expansion. 


assumption, which is reasonable for low values of Mos is a conservative 


one, i.e. Pr3/Pro is lower. 4 re - : 
Substituting: T,/T3 = [ + aa ei ] | + 5 Ms | 


| 
é 
and p) = p3, turns equation (4.6.1) to: 
jae) [i+ S (4.6.4) 


Solving equation (4.6.4) for My: 


BIG Be Ses ean at cat Set 


* « 
Bo cae ee re ee re ee 


- [1 +4a8-1 
31° VV co 


= ote 


M (4.6.5) | 


Where: 


NS 


DS nt ee Ry a ag eR a ne ee 


2 

y,c1 AoM 

a gee (1 toMy”) — (4.6.6) | 
3 


The subscript I shows that the computation was made from the inlet 


“~~ ot 


direction. 


By assuming again that po = P3> 


p+ Pr3/P 
gees (4.6.7) 


| 
! 
A4.6.2 Pressure 
| 
| 
: Bre Pra/Pe 


aS ee 


Pr3i = Pro Teva 8 (4.6.8) 


Where subscript I is as defined previously. 


48 


ia wa seh AR a aed ad Se 
z 2 . 


Fae Fer eee 


| 
| 


and therefore: | 


i 


Melee a nn 


oe 
Sees Ep BS 


ETRE N toe fete ent a ee 


ce TNS See ee Re Re on, 


A4.7 Location of Normal Shock Wave 

The main problem in computing the Mach numbers and the total pres- 
sures at the inlet arises from the fact that the exact location of the 
normal shock wave is not known, and must be found. The way of solving 
this problem is as follows: 

First, solve for two extreme conditions by assuming that the normal 
shock wave is located at the throat and at the exit of the inlet, respec- 
tively. After kaowing the lower and the upper values for Mars Pry: 
iteration can be made to find the exact location cf the normal shock 
wave, The criteria for this iteration is matching of values for Mas 
Pr3 from both the inlet and the nozzle directions, j.e.: 


(4.7.13 
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APPENDIX 8: TRAJECTORY EQUATIONS 


Bl. Atmospheric Functions 


Best fit curves were calculated for basic atmospheric functions, 


pressure, density, temoerature and viscosity of air as a function of 


the flight altitude. The basic formula which was used for this process 


is as follows: 


F = A exp(-B X 107° n°) (1.1) 


Where h is the altitude in meters, ana A,8,C are numerical parameters. 


The appropriate atmospheric functions are as follows: 
-6 nt 03) (4,2) 


Py = 1.03322 X 10° exp (-59.148 X 10 
) (1.3) 


Py = 1.228845 exp(-29.0144 x 10°S n ! 
Ty = 288.26 exp(-13.232 x 10°F n 1.0709) 


When: h = G ~ 17,000m 
Ty = 217.24°, when h = 31 ,000-32,000n {}.4a) 

Up = 1.793 X 107° exp(-45.1374 x 167 pO-59ES) 

When: ob s G-11,900m 

z (75a) 


Cc 
“ when A> LL,oon m 


Ug * 1.41724 % 10 

In these farmutze, the pressure (ey) has dimensisas of Kg/me, tne 
2 te z : : 

density (Gy) is given in kg/m”, the temperature rr is given in OV» 


and the viscosity (ug) is given in kg/(m-sec) for: A-sec/m°). 


B2. Drag 
Be.t Cowl Drag Coefricient 
2.1.1 It was found chet the cowl drag coefficient has a strong 


influence on the resuits. Therefore, a new model for this cow) drag 


coefficient was developed. The model was based on a theoretical 


development done wreviously by Prof. T. H. Gawain [9]. The main 
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difference between this development and the classical theory, is that 
the boundary conditions are applied at the body surface rather than 
along the axis. 

2.1.2 The model, which originally was developed for simple cases 
(cones, etc.) was modified to fit the shape of the projectile, illustrated 
in Figure 1.1. 

2.1.3. The modified program is listed in Appendix G. In the combined 
program (TRAJET) an interpulation procedure was used as a subroutine in 
order to simplify the calculation process. 


82.2 Base Drag 
; : After checking the influence of the nozzle exit area (Ag) on the 


performance, it was decided to aliow Ag to reach the maxim value 


possible (A.)s in order to reduce base drag. In this case, the base drag 


ae is negligible. 


tet ; 


i B2.3 Skin Drag Coefficient 
| be _ The Reynolds number is well known to te: 
|| Pqlot | 
i Re, = a (2.3.1) 
jae The transition Reynolds number of 


Re* = 2«10° 


is usually taken as criterion for transition between laminar flow 


eee ae 


(lower values) and turbulent flow (1. , ... values). The incompressible 
| 
i jaminar skin friction coefficient is related to the Reynolds number 


as follows: 


= 1.328//Re, (2.3.2) 


1 Rt LT NNT. | pcr veitniterme ne . 


Cyst 


Where the subscripts DS,L stand for skin drag coefficient, and laminar 


flow, respectively. On the other hand, the turbulent skin friction 


{ 
rY 
t . 
| coefficient (Cys q) is calculated indirectly from the formula: 
t 51 
t 


Wo ae fis DES es ee ee ee =~ oe See PS Ee A a, Sa IT dee YE ER erties cre hae 
at 
2 ; : 
L 
ek 
| i 
p5,7 10d, g (Cog 7p Rey) = 0.242 (2.3.3) | 
y 
The computation is done in subroutine CALDC, which works in a similar | 
way to subroutine CALCM which has been described earlier concerning the . 
1 | 
calculation of various Mach numbers. ‘ | 
a I 
B2.4 Wing and Fin Drag Coefficients ' he 
{ 
| The wings and the fins of the projectile, also contribute to | 
| drag. Basically, each of these drag coefficients contains two parts: | ; 
l 
| -Wing/fin wave drag. | r 
-Wing/fin friction drag. ie 
. ‘i { tea 
' al 
2.4.1 A psuedo 3-dimensional model was chosen to simulate the wave | . 
ni coefficient, The basic formulae used in this calculation are: | 
i ie } 1 
i) TS i 
y.11 ie -1 meas 
wotesess =I 3 24). “1 2. 2.4.1 Zz 
v(Mo) = Ya tan J Y,r1 (My 1) - tan J My ] ( ) ee 
| i 
{ i 
| fath 9 val (vat) 
Po/Pro = (1 +A My) 28 (2.4.2) 3 
| 
L' The effective span could be taken as (see fig. 82.1): 
bi=b-% (2.4.3) | | 
Substituting: tam = 2/c , simu = 1/M gives: | 
| : 
| bi sb - —S—— 2.4.4 3 
ME ran el 
Myo-1 i 
ie The drag coefficient would, therefore, be: 
i cw 2 fPor/Pro ~~ —Poz/Pro | tb: (2.4.5) 
- pee Me Pg/Pro Po/Prg |] AL 
= 
i The way ths simulation uses the above equations could easily be 
\ understood when looking at the formulae together with the flow chart of ; | 
the appropriate subroutine (Appendix C). Interference drag is ignored. | 
| | 
52 
ft | 


Pi 
7 a mye ner tte arya me Ne ean Re wee tans 


ie os ae hee Ree eee « ~ 
PO RRS a BE ole i St ik ah i nS Pe ee Laas - 


ne OE NO SET ECE URN eerste pte en : ee een) 


“ 


8 
ff 
: 
is o 
f, > 
it img 
A - 
x w 
“i ” 
. = » 
r os 
. at 
fl KH 
r oe = 
{ _ 
g oOo Hn 
I — 
iy oe a 
E (Sere | 
, -~ & 
4 ~~ Oo 
oc = 
e& v 
vw Bu 
—& vv 
iS) 
(Fa) Ld 
t ~4 
L<] 
a Oo 
“ww 
Nw 
am & 
y 
S 
PJ 
oo 
— 
pees 
' 
~ ee AR ee meer tem * pee * * (atpeinenetetihden beers ananereemnanen tetas i 
wa ee Seat ¢: a eae Stee 
a icine SiS SS ah Le fies a ae st as A Ee ats aaa ls Os ie cae ie ats eee a ae pe eee) OES eS Pr ee oe hs a en rr a | 
POR S : Ae? nd # 


ae ee 
CES el tap catean “ats Be 


2.4.2 The friction drag coefficient was calculated using the 


existing model for skin drag coefficient. (See Section B2.3) 


B2.5 Calculation of Drag 


Define the dynamic pressure (9) as follows: 
G's pQg? (2.5.1) 


When a is in units of kg/ (m.sec2), (or: nym?) . Also, define the following 


geometrical units: 


A_ = RS; S, = 2nkL (Ace) 


Where R,L are the radius and the length of the projectile, respectively. 


Similiarly, for the wings or the fins: 
Suny = nbc (2.5.3) 


tr 
Where n ts the total number of wings/fins (a value of 8 was taken fur a), 


and b,c are the span and the chord of the wing/fin (see fig. 82.1) 
Consequently, the drag (D) is given by: 


D=4q A, Con + 3p Cas * Su (Cong + Cour) (2.5.4) 


82.6 Orag Coefficient of a Conventional Projectile Without 
Propulsion 


The program has an option to calculate also the trajectory of a 


projectile without a pronulsion. 
The formulae used to clacuiate the drag coefficients in this case are: 


Nose Drag 
Cpy = (0.083 + 0.096/M,") (a/10)* 89 
is the cone half angle. 


2.6.1 


Where ao 
2.6.2 Base Orag 
Cop = (0.6837 - 0.3165 M + 0.0525 mM?) (Z/m) 
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The projectile is a conventional round. 
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2.6.3 Skin Drag 

Skin drag is calculated as discussed in section 82.3. 

é.5.4 Drag 

D = Drag =q A (Coy * Cpg) * So “ps 

63. Boosver 

The projectile has an initial muzzle velocity of 2500 ft/sec. Part 
of the combustor volume can be used as a booster to accelerate the 
projectile even more so that starting the ramjet will be easier. 
Define exhaust velocity (U,) as: 
(3.1) 


Ve i ‘50,8 9 


Where I is the specific impulse of the booster's fuel (in sec) and g 


$p.8 
is the acceleration of gravity (in m/sec}. 
From Newton's law (6, p. 323] 
. Z a du 
F = th, U, (m, fiat) ag (3.2) 
Where m_ and the are the mass of the projectile and the mass flow of 


P 
the booster respectively. 


it 
dU = th, U, -—e— (3.3) 
-Consequurntiy: 

[og = figt 

QU = U(t) - UO) = -U,tn —_——— (3.4) 
Where + is the booster burn time. Hence, 

m 
Aye U0 ——P 4.5 
e m, - Mp ( . ) 


AU is the change in initial velocity due to the booster, where Me is the 


mass of the booster. 
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B4. Oynamics 
The flat earth trajectory with drag and thrust is wel? known. 


The differential equations of motion are: 


42 


dy Se ps : 

re aq + (F-D) sind/m, (4.1) 
2 
d’x 

ie = (F-D) cose (4.2) 


where y is the altitude, t is the cime. g is the acceleration of 


gravity, F is the thrust, 0D is the drag, 9 is the clevation angie, 


6 satan etnah mall 4 


i mM, is the projectile mass. 
For numerical solution of equations (4.1) and (4.2) a finite difterence 


form can be used as follows: 


a 
| 
2 
| X49 > (F-B) cose At /m, * 2x54; - x, (4.3) | 
fe 
| : . Yjag = L-g + (F-D) sinosm,} at® + 2y5,) -¥, (4.4) | 
ba The values for j = 1 are from initial conditions, i.e. 
it 
5) ne xy = 0, yy] = 0 (4.5) 
|" ne 
: Xo = Up * cos@ » At, Yo * Ug * Sing» Ot (4,6) | 
, P = y t 
|: § = arctan wi SS] oe (4.7) 
: jt2 ~ *j+) 


Having calculated X42 and “440 from (eq. 4.3, and 4.4), one obtained new 


values for trajectory parameters using (eq. 4.7). Also thrust, drag 


and projectile mass are updated for j+2. 
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APPENDIX E: COMPUTER PROGRAM LIST OF SYMBOLS 


PROGRAM EQUATIONS UNITS MEANING 
SYMBOL SYMBOL 


GEOMETRICAL SYMBOLS 


Poveda OIE ere 


AREF A. mé Reference area 
| AJ A; me Area at station j 
' AIC Ae mé@ Area behind a conical shock wave \ 
i 
ie Adi* Ray mn? Area ahead of a normal shock wave 
| . “a AJd2* Axo me Area behind a normal shock wave i 
( | 
A30 Aa m* Initial area at station 3 Coa 
\ 
| ATAd AL/A, - Area ratios 
! | 
- | 
eae | 
i L3 ly m Length of combustion chamber : 
£ : Ls 
ae LPR by m Length of projectile 
| RPR Ry i Radius of projectile, ‘ 
ALFA & deg. Inlet cone half angle 
an 
ee 
Ne 
a | % 
' 
i! 
oa 
. te 
is 
43 
“ 
| ; 
, , nn 
eau 
Ve | “hen:  J=] 2.9, the shock wave is at station 1 2 real case, respectively. 
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| SYMBOL 


fy i 
z GAZ 
i : GAT2 
# GA2 
l i 
| : 
my 
a 
i 
1] 
ie 
j.! 
ad 
a 
oe 
Ly 
ime 
{| 
i 
HY, 
Vy, 
oo 
i! 
“| im , 
# fies ot 


EQUATIONS 
SYMBOL 


UNITS 


MEANING 


ATMOSPHERIC SYMBOLS 


% 
kg/m" 
kg/m? 
N.sec/m® 


m/sec 


7 


re rt te rents nrg or ar = te et ern nen cemeteries pre ee ts ee Et termes 


Static temperature at altitude Y3 


Static pressure at altitude ¥3 


Air density at altitude Y3 


Air viscosity at altitude Y3 


Projectile muzzle velocity 


Projectile initial mach number 


Gravity (9.807) 


Air heat capacities ratio (c,/cy) 


(yt! )/2 
(y,-1)/2 
(yg */[2(y,-1)] 
tal (Ya7)) 


107 


i ‘ 
i tile oe 
i 
| 
| 
4 

| 


q PR a ae epg é z 
( = = i aa a OT weber SL Mee OP RR em me le Boe: 
“4 if 
7 
i 
! i 
a PROGRAM EQUATIONS UNITS MEANING - 
M SYMBOL SYMBOL, 
| COMBUSTION CHAMBER'S SYMBOLS f 
be RHOF Oe kg/m? Fuel density 
H ' 
ia ETAT ny - Burning efficiency 
A,N A,N - Burning rate parameters ; 
| WA we kg/sec Air mass flow 
| WF We kg/Sec Fuel mass flow 
WT wr kg/sec Total mass flow 
| F F - We/ Ws 
rit ROOT ‘a m/sec Burning rate : 
of CSTAR cr m/sec VERGE / T(when:t = P(y,)). 2 
i f* 14 Fo | 
th 3 RF Re m/°K Hot gas constant] + R(J/mole/ ] 
- MW(kg/mole) +g (m/sec*) hows 
S GF \¢ 2 Hot gas heat canacities ratio (c5/c,,) 
Ty GF - : (ygtl)/2 l. 
e GF2 - - (vg 12 a: 
i GF}2 - - (vptT)/ 02 (r—-1)] a, yes 
a ‘ fw! 
. GF3 - : ¥¢/ (Y¢-1) Ds 
\ ! a - ; 
Par | 
io a 
. 
Kes i i 
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sae ; : 
i po 
i Be sire 
~~ | | 
‘ aan 
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PROGRAM 
SYMBOL 


TJ 

TT) 

PJ 

P TJ 

Md 

PTJT, PTI2 
MJT ,MJ2 


T3M 
M3N, M31 


EQUATIONS 


SYMBOL 


Te ee ae 
Fi TTP 


« 


UNITS MEANING 


THERMODYNAMIC SYMBOLS 


Se Static temperature at station j 
0% Total temperature at station j 
kg/m? Static pressure at station j 
kg/m? Total pressure at station j 
- Mach number at station j 
As above with Ad?, AJ2 
kg/m? Maximum 13 available 


M, calculated from nozzle and inlet 


direction, respectively. 
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PROGRAM EQUATIONS UNITS MEANING (Pd 
SYMBOI. SYMBOL 
LOSSES SYMBOLS 
! | PID} 1) - Boundary layer lcsses — 
| PID2 Ty" - Subsonic diffuser recovery | 
PIN a - Nozzle losses Be 
a0 ! 
4 - Te - Conical wave losses 
| | - Ts - Normal shock losses | 
| { 
- T. Lo Expansion losses | 
| 7 Th - Heat losses | 
re 
: | 
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PROGRAM EQUATIONS UNITS MEANING Lg 
SYMBOL SYMBOL. | 
RAMJET PERFORMANCE SYMBOLS | 
| Cr G i Thrust coefficient i 
2 cel 
sf Thrust F N Thrust ; 
4 {or kg) | 
: | 
' ISP i N/kg.sec Fuel specific impulse 
| P (or sec) | 
| Poy 
| i : 
e SFC SFC kg/hour/N Specific fuel consumption ce 
I I } 
[ = 
i I: 
fa ; 
a i 4 
i, 
cee va 
ce iy 
i . 
aa | 
Es ie 
7 Paes 
a ; ae 
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i om 
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PROGRAM EQUATIONS UNITS -- MEANING 
SYMBOL SYMBOL, | 
| 


TRAJECTORY SYMBOLS 


Drag D N Orag 
WPR Mp kg Mass of projectile 
WB We kg Mass of booster 
| DELU AU m/sec Change in initial velocity due to 
| booster 
U U m/sec Ug + Au 
s | DELT AT sec Change in time | 1 
| | ; TI t sac Time a os 
: TOB top sec Time of burning ! 
a : TETA @ deg. Gun elevation angle | , 
| ; TETP 85 dey. Projectile second cowl angle ; ce 
: DRAG COEFFICIENT, . 
cDN ea . Nose drag coefficient E* 
| CDWW coWW - Wing wave drag coefficient | | . 
y | CDMF Gaui - Wing friction drag coefficient {of 
i. es Che - Skin drag coefficient (Taminar/ turbulent) { I, 
|" CDB Cog - Base drag coefficient | : : | 
: | 
e oS 


— 


jaca, ee ee ha I a 


Jotceens 


OT sap ae a ae ~ os = P 
ih ; Sree corer ee te tre ee Se PD el mete ce on 


eet ae 


PROGRAM EQUATIONS UNITS MEANING 
SYMBOL SYMBOL 


MATHEMATICAL SYMBOLS 
v - 3.14159 


ms) 
rar 


IPR - - Printing parameter: > 0 combustion 
results together with trajectory 
(on different files): | 


= 0 working results only; 


=] also reasons for not running 


2 also full reasons for not running 


= 3 also Joop on mach number (CALCM) 


= -] trajectory prints only 


ITRA - - Loop paramerer: 


eae 


x. 


= +] single value for Al Ays Ag/A., 


= -] loop on Ap/A., Ac/A., and 


s 


print summary, only. 


IL00 - - Check parameter : 


< 1 regular run 


uf 


> 1 doesn't have a solution 


oh amet OS cr 


IL, IL0 - - Trajectory printing parameter 


* 


: (prints every ILO point}. 
Le IRAM - - Ramjet parameter: 


i = 0 ramjet in operation 


i = } projectile without propulsion 
i 


XMO Xu - Stopping mach number ) ; 
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‘ 
nS epricamermip: San water eae ee gare 


SO a 


PROGRAM 
SYMBOL 


INIT 
BURN 
TERMO 
INTER 


NOZZ 
CALCM 
RESUL 
CHOKE 
INLET 
CONE 
THROAT 
NSR 
DIFFUS 
MIN 


-CORVAL 
EXPAN 
HEAT 
CHECK 


TRAJ 
AT™ 


EQUATIONS 


SYMBOL 


a a ts Sl a i 


eA ee ee ee 


MEANING 


SUBROUTINES 

Computes initial conditions 
Computes combustian chamber's performance 
Thermodynamic comand subroutine 
Computes, by interpulation, thermodynaic 
conditions, or cowl drag coefficients. 
Computes nozzle performance 

Computes mach number indirectly 

Computes ramjet performance 

Checks if nozzle is choked 

Command subroutine for inlet 

Computes conical wave loss 

Computes boundary layer loss 

Computes normal shock loss 

Computes subsonic diffuser performance 
Command subroutine to compute situation 
when normal shuck is at point 2 

Command subroutine to compute situation when 
normal shock is at the correct place 
Computes Tosses due to expansion into the 
combustion chamber 

Computes heat losses at combustion chamber 
Check for pressure capability 

Conmand subroutine to compute trajectory 
Computes atmospheric conditions as a function 


of altitude (y) 
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PROGRAM EQUATIONS MEANING 

SYMAOL SYMBOL 

BOOS - Computes booster performance 

DRAGG - Computes drag 

CALCD - Computes skin drag coefficient, indirectly 
DRAGH - Cowl drag command subroutine 

DRAGHW - Computes wing/fin wave drag coefficient 


(command subroutine) 

PRANT - Computes Prandtl-Meyer angle (v) 
from a given mach number 

APRANT - Computes mach number from a given Prandtl- 
Meyer angle (v) 


PRES - Pressure ratios formula 

DRAGPR - Computes drag coefficients of a projectile 
without combustion 

DYNA - Computes the dynamics of the projectile 

DATTA - Initial data 

PRILOO - Prints ramjet performance 

PRIN - Prints detailed values when does not fine 
solution 
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APPENDIX F: COMPUTER PROGRAM USERS GUIDE 


Input Data 
AOAR, ASAR, A1AQ, A2A0, A3AR 


ALFA, TETA, TETP 
IPR, ITRA, XMO, TRAM 
PIDI, PID2, PIN 
Options 
ALFA = Inlet cone half angle 
TETA = Gun elevation angle 
TETP = Second cowl angle 
IRP = 0 clean print of RAM + TRAJ 
= 7, 2, 3 more details on RAM 
= «) TRAG only 
ITRA = 1 works on one set of data 
= =) loop on AORZAS 
TRAM = 0 Ramjet 
= | Projectile without rropulsion 
Execution Commands (For use with 18M 370) 
F2.1 Opening Consands 


LOXXXXP 
Password 


GLOBAL TXTLIB FORTMOG2 MODZEEH 


F2.2 Compilation 

FORTG] TRAJET 

F2.3 Run on Toning] 

FILEDEF 02 DISK TRI D(RECFM F BLOCK 132 PERM 
FZLEDEF 03 DISK DRG D(RECFM F BLOCK 132 PERM 
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Rat efter tree ei ae 2 BE RAO We RNS aye et se” 


ree Bae eae me ee + 


oP a 


eenevesrereminyeagpereper eine cnet: nee a 6 | tn An PN ert Se IAI AEG SNOT 


f 
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| 
t 
| 


7 on So a ee a aoe cae eas pe ee nee eee a ig ST FO ee eee ee ee is 
eH \ 
it 
; | FILEDEF OS OISK INP O 
' | FILEDEF 06 DISK CMB D 
EXEC RUN TRAJET 
XEDIT CMB OD | 


(or: 


(or: 


XEDIT TRI D) 
XEDIT DRG 0) 


' 
- remy 


Note: CMB D will cuntain the combustion process results. 
TRJ 0 will contain the trajectory part. 


ORG _D will contain drag coefficients. 


~~ 


Se ee ee 


F2.4 Hard Copy 


FILE 

PRINT CMB OD 
PRINT TRJ OD 
PRINT ORG O 
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G4. Results from AERO/COWL 


G4.1 The symbols used in programs AERO and COWL are defined in AERO 


and are presented in Figure G4.1. The values of a, bys Cys bs» Co are 


normalized with respect to n- NA, NB, NC, NBos Nl» are the appropriavs 
numbers of points used in the program. The cowl angles G, +a were 
selected as 20°, 9.5°, respectively. 

To create ficw at the first cowl in Figure G4.1 which ts the same as 
for a ramjet inlet, an extenston to the body was used. The extension con- H 
sists of the cone of length a and cylinder of length by in Figure G4.1. | 
Yhe cone angle, By> is 8 degrses. The value of b, was varied unvil the 


pressure at the first cowl was equal to ambient; a value of by equal to 


ry gives this condition. 
G4.2 The normalized values of the various variables which were 


selected are as follow: 
a= 7.it , by = 7 44, 20.2 4 dg * 2.13 , cy * 1.06 units : 
~ DL = 0.0736, ~ DL, = 0.177 units. 

Yhe appropriate numbers of noints are: 


NA = 203, NB = 200, NC = 6, WB, = 61, NC, = 30 


ot Sate ere one RR me ee 
? 


The dimens{onal values are: 


c, * 6.40, b, = 4.22, Cy * 2.1 inch 


3 


rm," 1.98, ry = 2.13, ra ™ 2.48 inch 
ec 


G4.3 It was found that the cowl drag ccerficient is sensitive to the 
magnitude of projected cow] area. For example, another combination of 


these vartables: 
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nt mentee tO TT ey OF Re mm ee %. 


ee 


aad Bae cats a 


ee 
7 _ epee ae emma z ” 
= eS ee a oe eee = = 
me er re ee a ae SP 


az/7.ll, Dy = 71, c * 0.38, be = 2.16, Cy * 1.46 


NA = 196, 8 = 193, NC = 11, NBp = 60, NC, = 40 
“DL = 0.138(20°), -DL, > 0.243 (9.5°) 


Gives higher values of the cowl drag coefficient: 


Mo 3.0 2.3 


Therefore, attention was made to select an inlet shape as smooth as 


G4.4 AERO can produce alse graphical results. A typical example is 


illustrated in Figure G4.2. 


CA.5 Computer Program Users Guide 
Use opening commends and comptlation similar to that defined in 


F2.1, F2.2. 


The routine to run AERO on the terminal is defined in pragrem 


itsel’. 

The routine to run COWL on the terminal fs as follows: 
FILEDEF 02 DISK CUT D (RECFM F BLOCK 132 PERM 
EXEC RUN COWL 
PRINT OUT DO 
XEDIT OUT O 


{ optional By 
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APPENDIX H: RESULTS 


Hi. SUMMARY 


Ay/A,. Ag/A. Ay /Bg Agi Mo Ay/A, ey 9 
0.25-0.40 0.25-0.40 0.47 0.887 0.426 98 15 
25 . 
35 | 
4c 
45 
65 
80 | 
G.25-0.40 0.25-0.40 0.50 0.827 0.826 9.5? a 
25 a 
45 oe 
65 
80 
0.25-0.40 0.25-0.40 0.42 0.827 0.426 9.5 45 mie 
6.58 
0.25-0.40 0.25-0.40 0.47 0.75 0.426 9.6 45 | 
0.827 tan 
0.91 ; 
0.25-0.40 0.25-0.40 0.47 0.827 0.27 9.5 45 a 
0.32 | 
0.47 
0.25-0.40 0.25-0.40 0.47 0.827 0.426 6.5 45 l 
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